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ABSTRACT 

The various four-fermion production channels in e + e~ annihilation are discussed and the CC11 
process e + e~ — > fx ft f% fit f^ e i ls studied in detail. The cross section d 2 a/dsids2, with s 1; s 2 
being the invariant masses squared of the two fermion pairs, may be expressed by six generic 
functions. All but one may be found in the literature. The cross section, including initial state 
radiation and the Coulomb correction, is discussed and compared with other calculations from 
low energies up to t/s = 2 TeV. 
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Figure 1: The double resonating CC3 contributions to off shell W pair production: crayfish 
and crab. They are part of all the processes of table 1. 

1 Introduction 

A measurement of W pair production around and above the production threshold will be one 
of the main tasks of LEP 2 as well as of a future high energy linear e + e~ collider [Q. 

Immediately after their creation, the W bosons decay and four-fermion (4/) production is 
observed: 

e+e" -> (W~W + ) -> 4/. (1.1) 

This double resonating process proceeds via the diagrams shown in figure |l[ In addition, there 
are a lot of Feynman diagrams with the same final state, but different intermediate states, 
which are single resonant or non-resonant and often called background diagrams. Their number 
and complexity vary in dependence on the composition of the final state. The total numbers 
of Feynman diagrams for the different channels are shown in table 1. 
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Table 1: Number of Feynman diagrams contributing to the production of two fermion doublets. 



One may distinguish three different event classes, all of them containing the CC3 process[) 
(i) The CC11 process. 

The two fermion pairs are different, the final state does not contain identical particles 
nor electrons or electron neutrinos (numbers in table 1 in boldface). The corresponding 
eleven diagrams are shown in figures 1 and 2. There are less diagrams if neutrinos are 
produced (CC9, CC10 processes). 

1 In ||, a slightly different classification has been introduced; the relation of both schemes is discussed 
in §. 
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Figure 2: Single resonant contributions to off shell W pair production: a w-deer and a d-deer. 





(ii) The CC20 process. 

The final state contains one e together with its neutrino (Roman numbers in table 1); 
compared to case (i), the additional diagrams have a t channel gauge boson exchange. 
For a purely leptonic final state, a CC18 process results. 

(iii) The mix43 and mix56 processes. 

Two mutually charge conjugated fermion pairs are produced (italic numbers in table 1). 
Differing from cases (i) and (ii), the diagrams may contain neutral boson exchanges^ 
(see also table 2 in [|| where 'neutral current' type final states are classified). There are 
less diagrams in the mix43 process if neutrinos are produced (mixl9 process). 

In this article, we will investigate the simplest topology, case (i), which proceeds via the 
diagrams of figures [I] and 

e + e" - Kftftfl, fi + /a, fi ± e. (1.2) 

As far as observable final states are concerned, the following reactions are covered by this 
study: 



(i) (lu/p) + (tu t ); 

(ii) (Wi) + (q u q d ), I 

(iii) (ud) + (cs). 



We parameterize the seven-dimensional four particle phase space as a sequence of two 
particle phase spaces: 



dT 
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7T JX{s,s 1 ,s 2 ) J\{sx,ml,ml) J\(s 2 ,ml,mj) 



256 



si 



S2 



dsi ds 2 dcosOdfli c?f2 2 , (1.3) 



2 We exclude the Higgs boson exchange diagrams from the discussion. In fact, the methods developed here 
are applicable also for the study of associated Higgs production || . 
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with the usual definition of the A function, 



X(a,b,c) = a 2 + b 2 + c 2 - 2ab - 2ac - 26c, (1.4) 
A ee X(s,sx,s 2 ). (1.5) 

In ( |1.3|) , the rotation angle around the beam axis has been integrated over already. Variables 
ki and k<i are the four-momenta of electron and positron and Pi,P2,P3,P4 are those of the final 
state particles ff , f% , f$ with pf = mf Q. The invariant masses squared of the fermion pairs 
are: 

s = (h + k 2 ) 2 , (1.6) 
si = (P1+P2) 2 , (1.7) 

S2 = (P3+P4?- (1.8) 

The W~ production angle 9 in the center of mass system is spanned by the vectors (p\ + P2) 

—* — — 

and fci. The spherical angle of p[ (ps) in the rest frame of the fermion pair [/1/1] ([/2/2D is 
flj (fi 2 ): dQi = dcos9id(j)i. In this frame, the p 2 (pi) points into the opposite direction with 
respect to p[ (ps) . The kinematical ranges of the integration variables are: 

(mi+m 2 ) 2 < si < - m 3 - m 4 ) 2 , (1.9) 

(m 3 + m 4 ) 2 < ^^(v^-v^T) 2 , (1.10) 

—1 < cos^, COS6 1 !, cos ^2 < I5 (1-11) 

< 01,02 <2tt. (1.12) 

We are interested in analytical formulae for distributions in invariant masses of fermion 
pairs. Thus, we have to integrate analytically over the five angular variables in (|1.3|) . The 
squared matrix elements have been derived with CompHEP and the angular integrations were 
performed with the aid of FORM || . 

The cross section contributions may be grouped into several classes of interferences: 

frcciiOO = / ds i ds 2 [^cc3 + ° 3f + ° vS + cr ff + a hh ] . (1.13) 

In ( |1.13| ), the CC3 process is accompanied by the interferences a 3 ? of deers with the crayfish 
diagrams, the interferences a v ^ of deers with the crab, the interferences of deers of one 
doublet, and the interferences a? 1 ? 2 of deers of different doublets. 

Each of these contributions may be described by the product of a coefficient function C, 
which is composed out of coupling constants and gauge boson propagators, and a kinematical 
function Q(s, si,S2)- The latter depends only on the three invariant masses and contains the 
dynamics of the corresponding hard scattering process. The Q functions, which are needed 
for a description of the CC11 process are shown in table f|. The double resonating CC3 process 
is known for long to be described by three kinematical functions 0: <?cc3> a = 33,3/,//. 
The double resonating neutral current NC2 process, e + e~ — ► (ZZ) — > 4/, proceeds via two 
diagrams of the crab type and is characterized by only one kinematical function {?nc20 103 ■ We 



3 Fermion masses are retained in the phase space definitions and in the basic kinematical relations. In all 
subsequent steps of the calculation, they will be neglected compared to invariants s, s%, and S2- 

4 Since not only Z pairs, but also Zj and 77 intermediate states exist one is faced in practice with an NC8 
process with the same topology like NC2. 
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also have to mention that the simplest final state configuration of a complete neutral current 
4/ process has 24 Feynman diagrams and requires, besides {? N c2, only one additional function 
£nc24 0- It will be shown below that a complete description of the CC11 process after the 
angular integrations is possible with adding only one further function Qccti- 



process 


C"CC3 


a 3 ' 




a ff 






£33 Q$f nff 

»CC3' ^CC3; ^CC3 


G 3f 

ycc3 


^CClli ^NC24 




ycCll; t/NC24 



Table 2: The complete set of kinematical functions describing the CC11 process. 



In the next section, notations and the formulae for the CC3 process are introduced. In 
section 3, the contributions from the background diagrams of figure 2 are presented. The 
photonic initial state and Coulomb corrections are added in section 4 and section 5 contains 
numerical results and a discussion. Several appendices are devoted to technical details of the 
calculation. 



2 The CC3 process 

2.1 The total cross section 

The double resonating cross section is: 



occsO; si, s 2 ) 



2 



7TS 



Cccs^ccsl 5 ; si, s 2 ) + CcL^ccsO; Si, s 2 ) + CSca</ccs(s; Si, s 2 



t nff 



The kinematical functions are known from 



£cc 3 (>; s i 5 s 2) 



^- [A + 12 (ssi + SiS 2 + s 2 s) 
48' 



[s — si — s 2 ) [A + 12s(ssi + s x s 2 + s 2 s) 



- 24 [ssi + s 1 s 2 + s 2 s] sis 2 £(s; s 1; s 2 ) 
1 



f2.r 



(2.2) 



(2.3) 



48 



[A + 12s(si + s 2 ) - 48s lS2 + 24(s - s x - s 2 )sis 2 £(s; Si, s 2 )] . (2.4) 



In the superscript, the '3' indicates to the matrix element with a three gauge boson vertex and 
'/'to the diagram with fermion propagator. The function C?cc3 nas been slightly simplified 
compared to @, 0. The logarithm in ( |2.3j ) and ( |2.4j) arises from the integration over the 
fermion propagator in the t channel (see appendix 0): 

rt \ 1 l s- si- s 2 + \/\ , , 

£(s; si, s 2 ) = -= In ? =. (2.5) 

V A s — si — s 2 — v A 
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Some properties of A and C are collected for the convenience of the reader in appendix [A]. 

For the coefficient functions, we choose generic definitions. At first glance this seems to 
lead to artificially complicated constructions. At a later stage, however, the usefulness will 
become quite evident. We begin with 

Cms = £ CS ee (E,s;F 1 ,s 1 ;F 2 ,s 2 ). (2.6) 

Vi,Vj,V k ,Vi=W 

The coefficient C l [ 22 nas been introduced for the neutral current NC2 process in eq. (11) of fBJ; 
with a slight change of notion for the current presentation, it becomes: 

Cf 2 l]ab (F 2 ,s 2 ;F l ,s 1 ;E,s) = ' 



(6^)2 D Vi (s)D v .( Sl )D* Vk ( S )D* Vl ( Sl ) 
L{fl V^LUl V k )L(fi, Vj)L(fl V5) 

+ R(f%, V^Rtfl V k )R{fl V 3 )R(f b 2 , V^N C (F 2 ) 

x [L(Fi, Vj)L{F\, V t ) + R(Fx, Vj)R{F Xi Vfi\ N^F,) 

x [L(E,V-)L(E,V k )+R(E,Vl)R(E,V k )]N e (E). (2.7) 

The labels E,F X , F 2 denote the corresponding members of weak isodoublets. For the charged 
current, all left- and right-handed fermion couplings are equal: 

L(e,W)=L(E,W)=L(F l ,W) = (2.8) 
R(e, W) = R(E, W) — R(F i7 W) = 0. (2.9) 

Below we will also need neutral current couplings where flavors will cause differences for the 
couplings. 

The boson propagator is: 

D v (s) = s-Ml + i^~sT v {s). (2.10) 
With these definitions and using the relation G^/V2 = g 2 /(8M^), one easily verifies that 

CSos = {G " M2w)2 P(si)p(s 2 ), (2.11) 



S\S 2 



where 



and 



" w - ;^4tSot x brw ' (2 - i2) 

r w w = E^#^- (2.13) 

j buy A 

The off shell width Tw(si) is used throughout this paper. It contains a sum over all open 
fermion decay channels / at energy ^Jsl, and BR(z) is the corresponding branching ratio. 
Equation ( |2.11| ) makes the presence of the two Breit-Wigner factors explicit. These are nor- 
malized such that limry^o p(si) = 5(s$ — x BR(i). 
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The other two coefficient functions are: 



C s 



nst 
U CC3 



C^ 3 (F 1 ,s 1 ;h,s h ;F,s f ) 



with 



K/,7) = e ^f, 

«(/,7) = 
#3(7) = gsw, 



^.z (6vr 2 ) 2 \D w { Sl )\ 2 \D w (s 2 )\ 2 D Vi (s)D* Vj (s) 
x 53(^)53(^)^(^1, H^)L 2 (F 2 , W)N C (F 1 )N C {F 2 ) 

x [L(e,Vi)^(e,^) + ^(e,V;)i2(e,^)]^c(^), (2.14) 

2 C£ 3 (F 1 ,s 1 ;e,s;F 2 ,s 2 ), (2.15) 

i= 7 ,Z;fc=W 

_^ Ke I 

(6vt2) 2 \D w ( Sl )\ 2 D* w (s 2 )D* Vj (s)D Vk (s f ) 

x L 2 (F 1? Vy)^ 3 (V;)L 2 (F 2 , W)L(/i, V fc )A^ c (F 1 )A^ c (F 2 ) 

x [L(e,V 3 )L(e,V k ) + R(e,V 3 )R(e,V k )]N c (E), (2.16) 



and e = gs^/ = V^tto, Q e = — 1, If = —\. 



gc w , 



w I ■> 



(2.17) 



2.2 The angular distribution 

It is quite useful to know not only the invariant mass distributions and the total cross section, 
but in addition also the distribution in the production angle of one of the W bosons: 



with 



and 



da CC3 
dcosO 



TTS^ 



CcC 3 £cC3(s;Sl>S2;COS#) 



+ Ccc 3 £cc 3 (s; Sl ' S2 > cos ^ + C cc3^cc 3 (s; si, s 2 ; cos6) 



t r.ff 



^CC 3 ( S ; S l, S 2;COS^) 

QcL(s;s 1 ,s 2 ;cos9) 
0cra(s;si,s 2 ;cos0) 



[ACi + 12s lS2 C 2 ] , 

s - si - s 2 )Ci - 
4sis 2 C 2 



4sis 2 [s(s l + s 2 ) - C 2 ] 



Ci + 



t?. 



c 2 



2s(si + s 2 ) + C 2 , 

- Si - s 2 - t v ) - S 1 S 2 . 

The scattering angle is contained in the denominator of the t channel propagator: 

t v = - (s — S\ — s 2 — y/X cos . 



(2.18) 

(2.19) 
(2.20) 

(2.21) 



(2.22) 
(2.23) 



(2.24) 
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With the aid of appendix |Cj, it is trivial to integrate over cos 9: 

1 r 1 

£cc 3 0, s i, s 2) = ^J^ 00 * 9 ^cc 3 ( s > s i ; s 2;cos^), a = 33, 3/,//. 



(2.25) 



The distribution (|2.18 ) is used for the description of W pair production in PYTHIA [[Tl|]. It 
also may be used for the calculation of moments like 



dsi ds2 d cos 9 cos" 9 



da 



dsids2d cos 9 



[2.26) 



as a check on the accuracy of Monte Carlo integrations [12| . The angular distribution is more 



interesting in the context of anomalous couplings of gauge bosons [13 



3 Background contributions 

We now come to the contributions from the eight background diagrams of figure || and from 
their interferences with the double resonating diagrams of figure |I[ For this purpose, we use 
a classification, which was introduced for neutral current processes in ||. A single resonant 
diagram with a virtual fermion / is called an /-deer. Strictly speaking, the deers are double 
resonating diagrams as crabs are: besides the one W resonance, they contain the s channel Z 
(or 7) propagator. One of the invariant masses, in which the diagram may become resonating, 
is 'eaten' by the fermion line. In the crab it is s, while in a deer either s\ or S2- In this 
language, the crab is an e-deer. In appendix |C| it is made plausible that in fact this observation 
may be used for a treatment of all the contributions on an equal footing. 

3.1 The crayfish-deer interferences cr 3 ^ 

The interferences of the crayfish diagram with the four types of deers are similar to the 
crayfish- crab interference: 



a 



3/ 



7TS Z 



n=l,2 



(3.1) 



The coefficient function is 



QlCll — ^223(^3— n> s 3-ni fni s n] e , s) 

Vj,Vk=l,Z 



(3.2) 



and (?cc3 is defined in ( |2.3|) and C223 in ( p.!6| ). 
3.2 The deer interferences 

The contribution from the square of two deers, which belong to one doublet F = (f u , f d ) 
occurs twice: 



a 



ff 



Vx v 
2 2^ 



TCS 



71=1,2 



fu fu fd fd 

CJ n J n 1 n J n 
CC11 1 ^CCll 



s) , (3.3) 



,/«/^ r u,d 
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with the coefficient function 



^CClT — ^422 (-^ri) s raj F m , S m ] E,s). (3.4) 

The C % l 22 ah is defined in ( |2.7|) . 

The kinematical functions are different for the pure squares of diagrams, where ( |2.4j ) is to 
be used, and for their interference. The latter may be expressed by known functions. Let us 
refer for a moment to the neutral current case. There, all couplings are equal and one gets 
the following relation: £ NC2 = Q uu + Q dd + Q ud . With Q uu = Q dd = g££ 3 , and Q ud = 0g? u , one 
gets the requested identity, 

^CClll 5 ! S l) S 2) — QnC2{s] Si, S 2 ) — 2 QcC3{ s ] s l) s 2), (3-5) 

with the neutral current function Q NC2 from JTOj] ; see appendix [D]. 



3.3 The deer interferences a^ 1 ^ 2 

There are four interferences among deer diagrams belonging to different doublets: 



TIS 2 



C^+Cf^ d ]g-i d 1 d (s;s u s 2 )+ \C™ + C™ G U ccUs;s u s 2 )\. (3.6) 



The coefficient function may be traced back to that introduced in ||: 

CM = £ C%f(E,s;f?, Sl ;f h 2 ,s 2 ). (3.7) 

The function is defined as follows: 
c zaa( E , s ', F u8i ;F 2 ,s 2 ) = jfte 



(6^)2 L> yi ( s )^( s )D y! ( Sl )D^,( S2 ) 
x [L(£, ^)L(£, Vfc) + R(E, Vi)R(E, V k )\ N C (E) 

x [L(Fi,Vi)L(Fi,Vj)L(Fi,Vk) — R(Fi,Vi)R(Fi,Vj)R(Fi,Vk)] N c (Fi) 
x [L(F 2 , V t )L(F 2 , V,)L{F 2 , V±) - R(F 2 , Vi)R(F 2 , Vj)R(F 2 , Vt)} N C {F 2 ). 

(3.8) 

Wherever a neutral gauge boson couples, the corresponding fermion e or f 2 etc. will replace 
the doublets E, F 2 in the arguments of the couplings L, R. 

For the kinematical functions, again a relation to a neutral current function may be es- 
tablished. The function £ N c24, explicitly given in appendix |D|, describes a sum of interfer- 
ences analogue to those considered here. Since the couplings are equal in that case, we 

have: £ N C24 = Q UlU ' 2 + Q Uld2 + g d ^ _|_ g d ^ _ Furthermore, using the symmetry relations 
gmui = gd ld2 = g uu,dd ^ guid2 = gdlU2 = g u,d^ wg m&y write g ^ = 2g uu,dd + 2g u,d^ ^ 

1 1 • n j j j • s-yuu,dd s->u,d • i 

relation allows to determine y C cii once fc/ C cii J s known: 

s-iuu,dd 1 s-iU,d /q n\ 
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By an explicit calculation, we obtained: 



3^3 



0ccii(s; s u 82) = -120s 4 ^^£(s 2 ; s, si)£(si; s 2 , s) 



— s 



1 + K^) +20s 2^2 Qn . 3 '"3* 



A 2 



30s d sis 2 - 



A 3 



1 



S — cr SiS 2 9 S + 0" 

— h lO 5 ^ 30s sis 2 ^ 3 



Si£(s 2 ;s,si) + s 2 £(s 1 ;s 2 ,s) 
s i^( s 2; s, Si) — S2-C(si; s 2 , s) 



12 



2 2\ 

s — cr) 



sa oSiSo 
1 + 12 T - 60s 2 -^ 



A 



A 2 



A 



where a = s\ + s 2 . 

For numerical applications, it is helpful to know the limit A — > 0: 



lim^ccii(s;si;s 2 ) 

A — >U 



64 



A 



160 



where A = si — s 2 . 



(3.10) 



1 + 9— +0(A 2 ), (3.11) 



3.4 The crab-deer interferences o v * 

There are four crab-deer interferences: 

° = 2 

™ n=l,2 

The coefficient functions are 

L-CCll 



^ [OjCll ^CCll( s 3-n; S, S n ) + Cecil Qccil ( s 3-n] S , S n ) . (3.12 



^ ^233 (-^3-n> s 3-n] E, S, f n , S n ), 



and the is defined in (|3.8|). 



(3.13) 



The kinematical functions are introduced in section B.3I. One should be aware of a factor i 



resulting in the limit of neutral couplings. Differing from the neutral current case, there is no 
charged current crossed crab diagram, which would give the same contribution as the crab 
diagram. 



4 QED corrections 

A complete treatment of QED corrections for four-fermion production is part of electroweak 



corrections and a quite ambitious task |14j. For the time being, one may try to restrict to 
initial state radiation (ISR) and, near the W pair production threshold, corrections from the 
Coulomb singularity. Both corrections influence the cross section significantly. 
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4.1 Flux function approach 

The cross section with QED corrections becomes in the flux function approach: 



d(TQED(s) 

dsids2 



s 



G(s'/s) <7ccii(s',si,s a ) + C(s')a C c3(s',s 1 ,s 2 ) 



non—univ / J 



QED 



S ,S!,S 2 



with 



and 



G(s'/s) = p e 



1 + S) + H{s',s) 



2a 



71 



In 



mi 



The virtual and soft corrections S and hard corrections H are [JL5 

r 2 



s 

5 (2) 

H(s',s) 



3 n a (tt' z 1 
4 A + 7 T"2. 



x 



IZERO + S {2) , 



a 



s 2 L 2 e + s±L e + s + S'i 



(2) 



pairs 5 



K 1 + 7) A + flW ' 



2 

tt^ 2 

7T 



h 2 Ll + h 1 L e + h +H. 



r(2) 

L pairs ' 



(4.1) 

(4.2) 

(4.3) 
(4.4) 

(4.5) 
(4.6) 
(4.7) 
(4.8) 



(2) (2) 

where the corrections SjJ ira and H^J irs due to fermion pair emission from the initial state have 



been determined in 16 



In the leading logarithmic approximation (LLA), the ISR is well defined. The complete 
0(a) corrections include additional contributions in case of t channel exchanges. They have 
been defined and determined for the double resonating CC3 (flTl) and NC2 ((T§|) processes. 
For a complete treatment of the CC11 process, this is not sufficient and has to be accomplished 
by a proper treatment of the crayfish-deer interferences; these are, however, small and the 
non-leading QED corrections to them the more. 

The Coulomb corrections apply to the double resonating diagrams 



C(s) 

(3 



a(2M ] 



w , 



13 



3m In 



Pw + 5-P' 



m 



W 



Mir - tM W T W , 



(4.9) 

(4.10) 

(4.11) 
(4.12) 
(4.13) 
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Besides the cross section, there is also experimental interest to know the radiative energy 
loss due to ISR. This quantity cannot be calculated in the flux function approach. However, 
the calculation of a similar quantity, namely the invariant mass loss (of final state fermions) 
is possible: 



(m ) = 1 [dcjcj I dS> — (l S '] d ° 

a J J s 2 \ s J dsids2ds' ' 



(4.14) 



where da / dsids 2 ds' is the contents of the curly brackets in ( fO 



4.2 Structure function approach 

Alternatively, the structure function approach may be realized with the following formula: 

i i < 

= J dx\ J dx 2 D(xi, s)D(x 2 , s)< cr ccll (xix 2 s, si, s 2 ) 



da. 



QED\b 

ds\ds 2 



+ C(xiX 2 s)a C c 3 {x 1 x 2 s, sx, s 2 ) >■ 



(4.15) 



The lower integration boundaries are: 



.) '-) 



S 
XiS 



To a very good accuracy, the structure function may be written as follows PD 



D(x,s) = {l-x)^ 2 - 1 ^{l + S)+H{x,s), 



and 



S 



H(x,s) 



e (3/4- 7B )/3 e /2 

r(i + f) 

-1(1 + x)^ 
2 V ' 2 



- 1, 



(4.16) 
(4.17) 



(4.18) 



(4.19) 



Ve 



< 

2a 

71 



1 



-4(1 + x) ln(l - x) + 3(1 + x) lnx - 4 



In a; 
1 — x 



— 5 — x 



(L e - IZETTA) . 



(4.20) 
(4.21) 



Within a given order of LLA, the structure function and the flux function approaches are 
related: 

G(s'/s) = J dxidx 2 D(xi, s)D(x 2 , s) 5 {^xix 2 - — ^ . (4.22) 
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The radiative energy loss is determined as follows: 

(E 7 ) = — [ ds 1 ds 2 [ dx 1 dx 2 D(x 1 ,s)D(x 2 ,s)^-[(l - xi) + (1 - x 2 )] — 
a J J 2 ds\ 



da 



ds 2 dxidx 2 ' 

(4.23) 



and the invariant mass loss (of final state fermions) is: 



(ra 7 ) 



— J dsids 2 J dxidx 2 D(xx,s)D(x 2: s)^-(l — 



da 



ds\ds 2 dx\dx 2 



(4.24) 



where da j ' ds\ds 2 dx\dx 2 is the content of the curly bracket in fl4.15|) . 

In table |3], we have collected the flags, which are used in GENTLE in order to define the 
details of the calculations. 



Flag 


FF 


SF 






ICONVL 





1 






ICOLMB 


0-5 




3 





IQEDHS 


0-4 




3 


2 


I ZERO 


0, 1 




1 





IZETTA 




0, 1 


1 


1 


ITVIRT 


0, 1 




1 





ITBREM 


0, 1 




1 






Table 3: Some of the flag settings in GENTLE for the flux function (FF) and the structure 
function (SF) approaches. The third column of flags contains recommended flag values for best 
estimates and the last one those chosen for the comparisons in table 6. 



5 Numerical results and discussion 



The numerical results are obtained with the Fortran program GENTLE [^TJ. The input quantities 
are chosen as follows: 

= 1.16639 x 10~ 5 GeV 2 , 
a(2M w ) = (128.07)" 1 , 
M z = 91.1888 GeV, 

M w = 80.23 GeV. (5.1) 
Further, for the QED corrections, 

a = (137.03599)- 1 . (5.2) 
The s dependent widths are Ty(s) = (y/s/M v )T V - For the Z width we use the value 

T z = 2.4974 GeV, (5.3) 
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thus circumventing a dependence of all the numbers on the top quark mass from a prediction 
of Tz in the Standard Model p2[ . For the W width the radiative corrections are small |23] 
and the following approximation is used: 



w = 9x ^£_w. (5.4) 



Finally, we need the weak mixing angle sin 2 9w- Again, we deviate from the pure Standard 
Model relation and use instead an effective weak mixing angle: 



na(2M w ) 



In table [|, cross section predictions are shown for the different channels of table |] over a 
large range of y/s. They are normalized by taking off the branching factors. The differences 
between the channels are seen to be minor above \fs ~ 160 GeV, where the double resonating 
CC3 diagram are dominating. Nevertheless, the predictions for the CC3 process and any of the 
CC11 processes start to deviate substantially above yfs = 200 GeV. 



d\Ui 


d 2 u 2 


BR(l)xBR(2) 


Iv 


Iv 


1/81 


Iv 


q d q u 


1/27 


q d q u 


Iv 


1/27 


qdqu 


q d q u 


1/9 


all 


all 


1 



Table 4: The branching ratios of the different types of final states of the CC11 process. 





CC3 


Lvlv' 


Ivqq' 


QQ'qq' 


ecu 


30 


1.4519xl0~ 7 


5.9295xl0- 6 


7.2478xl0~ 6 


5.0897xl0- 6 


6.1422xl0- 6 


60 


1.9358xl0~ 5 


4.0025xl0~ 5 


4.6879xl0" 5 


3.5441 xlO" 5 


4.1034xl0- 5 


91.189 


0.11225 


0.021329 


0.024551 


0.018975 


0.021715 


176 


16.225 


16.242 


16.243 


16.243 


16.243 


200 


18.578 


18.586 


18.588 


18.588 


18.588 


500 


7.3731 


7.3301 


7.3318 


7.3334 


7.3323 


1000 


2.9888 


2.9342 


2.9344 


2.9348 


2.9347 


2000 


1.5702 


1.5020 


1.5018 


1.5016 


1.5017 



Table 5: Total cross sections in nbarn for the different Born 4f production channels as func- 
tions of ^Js (in GeV). For this comparison, the branching ratios of table 4 are not taken into 
account in the single mode channels. 
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Tables ^ and [7] compare for LEP 2 and table |8] for linear collider (LC) energies the numer- 



ical predictions of GENTLE with those of the Monte Carlo programs WPHACT [|J , WWGENPV p5| 
and WTO the latter performs the numerical integrations with a deterministic approach. 
GENTLE may be used with both the flux function (FF) and the structure function (SF) ap- 
proaches, while the other programs are based on the structure function approach. Shown are 
CBom, &qed, (E-y) , and (m 7 ) as being introduced in section |. For the observables with QED 
corrections, the first two lines at each energy show that both treatments of QED corrections 
agree well. Further, the predictions of the different programs are in perfect agreement with 
each other. Having in mind that some numbers arise from ninefold numerical integrations, 
this is a nontrivial result proving the high level of technical precision reached by now with 
quite different techniques. The dependence of cross sections on the channels has with QED 
corrections the same tendencies as in the Born case. As an example for this we mention that 
(rrij) in table |7] shows, within the digits shown, no channel dependence at all. 







176 


190 


205 


O Born 


GENTLE 


1.8048 


2.0405 


2.0573 




WPHACT 


1.8048(1) 


2.0405(1) 


2.0574(1) 




WWGENPV 


1.8048(1) 


2.0405(1) 


2.0574(1) 


°~QED 


GENTLE, FF 


1.5067 


1.8148 


1.9022 




GENTLE, SF 


1.5046 


1.8124 


1.8998 




WPHACT 


1.5046(1) 


1.8124(1) 


1.8999(1) 




WTO 


1.5046(1) 


1.8124(3) 


1.8996(4) 




WWGENPV 


1.5045(1) 


1.8124(1) 


1.9000(1) 



Table 6: Results from different programs for the total CC11 cross sections e + e~ — > Af and 
e + e~ —>■ Af + 7 (both in nbarn) at LEP 2 energies for the QQ' qq' channel. The MC error 
estimates are indicated. 







176 


190 


205 


(E 7 ) 


GENTLE, SF 


1.1151 


2.1400 


3.2025 


WPHACT 


1.1149(2) 


2.1402(4) 


3.2030(6) 




WTO 


1.1149(3) 


2.1404(5) 


3.2023(8) 




WWGENPV 


1.1147(4) 


2.1400(5) 


3.2028(10) 


(m 7 ) 


GENTLE, FF 


1.1129 


2.1310 


3.1816 


WWGENPV 


1.1124(3) 


2.1316(5) 


3.1861(11) 



Table 7: Results from different programs for the average energy loss (E 7 ) and the invariant 
mass loss (m 7 ) (both in GeV) due to QED radiative corrections at LEP 2 energies in the Iv qq' 
channel. 
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500 


1000 


2000 


0~Born 


GENTLE 

T TT T /I T — 1 H T T 7 

WWGENPV 


0.81482 
0.81480(6) 


0.32609 
0.32602(6) 


0.16684 
0.16682(7) 


°~QED 


GENTLE, SF 
WPHACT 

WTO 
WWGENPV 


0.86950 
0.86956(9) 
0.86960(25) 
0.86956(14) 


0.36514 
0.36515(5) 

0.36530(35) 


0.18247 
0.18250(4) 

0.18247(13) 



Table 8: Results from different programs for the total CC11 cross sections e + e~ — *> Af and 
e + e~ — > 4/ + 7 (both in nbarn) at Linear Collider energies for the QQ' qq' channel. The MC 
error estimates are indicated. 



Figure |3] shows, over a large energy range, the CC3 process without and with QED correc- 
tions. The latter smear out the W pair production peak and produce a radiative tail at higher 
energies, crossing over the Born curve at ^/s ~ 240 GeV. The Z resonance in the crayfish 
diagram leads to an interference pattern around yfs ~ 91 GeV, with a set-on of the Z radiative 
tail until the W pair production starts to over-compensate it and to completely dominate the 
process. 

The same, but now for the CC11 process, is shown in figure (|. As to be expected, the 
differences to the foregoing figure are not pronounced, although visible. Both the Z and W 
pair peaks (seen in the set-in as minima in gqed/ &Born) are nearly untouched, while the tail 
of the Z between the two minima gets suppressed due to the background. 

Direct comparisons of the CC3 and CC11 processes are shown in figure [| without and in 
figure H with QED corrections. Quantitatively, the figures again do not differ too much. 
At the W pair production peak, the background influence is smallest, tending to negative 
contributions of few percents at large energies. Below the W pair production threshold, the 
CC3 process becomes heavily suppressed and background becomes equally important if not 
dominant. 



5.1 Conclusions 

We determined the so-called background contributions to the off shell W pair production 
cross section for the simplest final state configuration. Exploiting the inherent symmetries, 
we arrived at a remarkably compact presentation for the double differential distribution in 
the invariant masses of the W bosons. With no doubt, the angular distribution in the W 
production angle is also quite compact, as our expressions for the CC3 case indicate. Other 
angular distributions might be accessible as well; for the neutral current case, there are en- 
couraging indications for that [ 2"7| . Staying more differential could be welcome to colleagues 



who search for anomalous gauge boson couplings, while, on the other hand, the advantages 
of the semi-analytical approach are partly lost. Anyhow, we see no principal problems arising 
when anomalous couplings are included. A step in this direction has been done for the CC3 
process [p8H - 

Finally, we would like to comment on the prospects of applications of our techniques to 
other four-fermion final state topologies. For the production of four untagged hadronic jets, 
one has to add incoherently the qqgg final state, which should be straightforward. More 
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involved is the treatment of production of identical particles and/or of processes with gauge 
boson exchanges in the t channel. The first problem is due to permutations of final fermions in 
some of the diagrams, the second one by angular dependent gauge boson propagators, which 
occurs in e.g. the CC20 process and is related to gauge invariance violation . 

Time will show to what extent the technical obstacles mentioned will be overcome in the 
semi-analytical approach advocated here. 
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Figure 3: The CC3 process without and with QED corrections. 
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Figure 4: The CC11 process without and with QED corrections. 
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Figure 5: The ratio o/CCll and CC3 processes without QED corrections. 
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Figure 6: The ratio o/CCll and CC3 processes with QED corrections. 
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A Some properties of the auxiliary functions A(s,Si,S2) 
and £(s; Si, s 2/ 



Some of the properties of the two auxiliary functions A and £ are collected here. We begin 
with the A function: 



A = \{s,s 1 ,s 2 ) = X{s 1 ,s,s 2 ) = X{s 2 ,s 1 ,s) 
= s 2 + sf + si — 2ssi — 2sis 2 — 2s 2 s 



(A.l) 

The A function is symmetric in its arguments. From (|A.1|) , the zeroes may be read off easily. 
One is the virtual threshold (and integration boundary), 



the other one is unphysical. On the mass shell, the following limit holds: 



(A.2) 



sJX(s,M^,M^) = s/3(s,M$y), 
^M 2 W ) = f^k. 

The logarithm C is defined as follows: 

£(s;s 1 ,s 2 ) = C(s;s 2 ,s 1 ) 

1 s - si — s 2 + V\ 
\f~X s — S\ — s 2 — y/X 



(A.3) 



(A.4) 



This function is symmetric in the last two arguments. For vanishingly small invariant masses, 
it is: 



lim C(s; si, s 2 ) 



SI ,S2— >0 



S SiS 2 



(A.5) 



The other important limit is approached when A vanishes, which corresponds to the virtual 
threshold kinematics: 



lim £(s: si, s 2 ) 



A^O ' S — S\ — S 2 

In the on mass shell limit, the following simplification holds: 



+ . . . 



C(s;M^,M^ 



s(5 1-/3 



(A.6) 



(A.7) 



In the ultra-relativistic limit, i.e. when Myy becomes vanishingly small compared to s, one 
better uses 



lim as\M 2 w ,M 2 W) 



V S 



(A. 
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B Four-momenta and phase space variables 

In this appendix, we express the four-momenta fcj and pj in terms of the variables, which 
are used for the parameterization of the phase space: si,s 2 ,0,0i,(f>i. Then, the phase space 
integration over the squared matrix elements, which are expressions in the scalar products 
kikj, kiPi,pip n , i,j — 1, 2, I, n — 1, . . . 4, may be easily performed. 

In the center of mass frame, the initial state vectors k\ and k 2 are given by 

ki — (k , — k cos 9, 0, k cos 9~) , (B.l) 
k 2 = (k , kcos9, 0, -kcos9), (B.2) 



with 



k i 



k = ^, (B.3) 
J\{s,ml,rril) 

sk = ^ ' (R4) 



This allows to calculate k±k 2 . 

The products p\p 2 and p^p^ may be determined in the rest systems of the corresponding 
fermion pairs. For the 'compound' V\ (corresponding to the W~ boson in case of resonant 
production), the rest system R is defined by the condition 

P12 = P1+P2, (B.5) 

pf 2 = 0. (B.6) 

It is 

„ gi + mj - m\ 

Pi,o " 2^~' (B-7) 

o si - mj + mj 

The 3-momenta of the two fermions are pointing into opposite directions and we have: 



i-fli i-fli R ^JX(si,mlm 2 2 ) 

The components of the fermion momenta p\ and p 2 in the rest system of the fermion pair are 
(pf = Ko'^Ap&}' i=1 > 2 ) : 

pf — {pfo'Pil s i n ^1 cos 01)^2 s i n ^1 s i n 0i) P12 cos ^i} ) (B.10) 
P2 = \P2fi1 ~ P12 sm $i cos 0i, — Pi 2 sin#i sin 1; — pf 2 cos#ij . (B.ll) 

These relations allow the calculation of the product p±p 2 in terms of the integration variables; 
similar relations hold for the other fermion pair and the product P3P4. 

The remaining Lorentz invariant scalar products kiPj,pip^, etc., may be calculated by 
boosting all the final state momenta into the center of mass system. We will assume that 
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the 'compound' V\ moves along the positive z axis. Then, it is in the center of mass system 

(^34 = P3 +Pa)- 



S + Si — S 2 /_ , 

Pi2,o = 2v/i . (B- 12 ) 

P34,o = 2yi , (B.13) 

\M = \M = (B.14) 

With the aid of the above definitions the Lorentz transformation from the center of mass 
system (and back) may be expressed. The components of p\ transform from the rest system 
of V\ into the center of mass system according to: 

Pi,o = 7? 2 pfo + 7i2P^, (B.15) 

Pi,x = P?*, (B.16) 

Pi,, = P* v , (B.17) 

PM = 7i2P?o + 7i°2Pm, (B.18) 



where 



t „ 2 s f^ = £ + f ^ £l _ 



2 v /ss7 



|Pl2 



7i2 = ^ = 7^7=1 (B.20) 

and analogously for p 2 and for the other fermion pair. One should have in mind that for the 
'compound' V 2 the sign of the transformation of the z component is negative. 

As a result of all these relations, we may finally express the components of the four final 
state momenta in the center of mass system, p i = {Pifl,Pi, x ,Pi, y ,Pi,z}, i = 1, ... 4: 

Pi ={7i° 2 ^o + 7i 2 pf 2 cos^, pSsmfccos^, pf 2 sin 9, sin fc, 7M cos + 712^0} , 

P 2 = {7? 2 pJo - 7i 2 pf 2 cos ls -pf 2 sin 6 1 cos 0i, -pf 2 sin 6 1 sin l5 -7? 2 pf 2 cos 6 1 + 712^0} , 

Ps = {734^0 ~ 734P3I cos 2 , p& sin 2 cos 2 , pf 4 sin 2 sin 2 , 7° 4 pf 4 cos 2 - 734pJo} > 

P4 = {734^0 + 734P34 cos ^2, -P34 sin #2 cos 2 , -pf 4 sin 9 2 sin 2 , -734P34 cos 6 2 - l3APf, } ■ 

(B.21) 

With these relations, any of the scalar products may be expressed in terms of integration 
variables. It may be seen explicitly that in the center of mass system the (pi + p 2 ) and 
(P3 + Pi)i as we h as (ki + k 2 ), are independent of all angular variables while (pi — p 2 ) and 
(p 3 — p 4 ), as well as {k-y — k 2 ), depend each on only one angle 9i, 9 2 , 9, respectively. 
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C Phase space integrals 



The gauge boson propagators depend on one of the invariant masses squared s,Si,s 2 and 
are left for numerical integrations. Thus, the only angular integrals will arise from fermion 
propagators. In this respect, the CC3 process is distinguished by its simplicity since the final 
state integrations factorize completely and may be performed by tensor integration over 9i, <f>i 
and then the neutrino propagator in the crab diagram contains the only nontrivial angular 
dependence: 



D,. 



[(ki - k 2 ) - (pi + p 2 ) + (p 3 + Pa)\ 



si - s 2 



V\cos6 



1 



(C.l) 



Here, only the difference {k\ — k 2 ) depends on an angle (see the final remarks in appendix 
The one dimensional integration is: 



[A) t 



- / dcos9 A. 
2 



(C.2) 



-i 



Besides [cos n 6]q 



[1, 0, 1/3) for n = 0, 1, 2, the following integrals are used: 

C(s;si,s 2 ), 
1 



1 



S\S 2 



(C.3) 
(C.4) 



The background contributions contain either one or two (equal or different) fermion prop- 
agators. An example is: 



D 



h 



[{h + k 2 ) - (px - p 2 ) + (pa + Pa)Y 



si - s 2 



s + a/X cos 6*i 



1 



(C.5) 



which occurs in the second of the Feynman diagrams of figure |2|. Here, only (pi —p 2 ) depends 
on an angle. The other propagators are expressed similarly. Evidently, the angular integra- 
tions over 0i,02 are not influenced and those over 6,6 1 ,6 2 factorize and may be performed 
independently using the above table of integrals, thereby permuting s,Si,s 2 in the answers 
as needed. After a careful book-keeping, the background contributions are not much more 
involved than the CC3 case. 



D The neutral current kinematical functions 



The off shell Z pair production may be described by one kinematical function |jlO 

N 2 



£nc2(s; si, s 2 ) 



+ (3i + S 2 ) 

s — Si - s 2 



-C{s,s x ,s 2 ) - 2. 



(D.1) 



The process NC24, which proceeds via off shell gauge boson pairs, but also via single 
resonant diagrams of the deer type, may be described by two functions (? NC 2 and (? N c24- The 
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latter is 0: 



C?nc24(s; si, s 2 ) 



SS1S2 x — <j £(s 2 ; s, si)£(si; s 2 , s)4s 



ss 1 (s - si) 2 + ss 2 (s - s 2 ) 2 



+ (s + Si + s 2 ) 



+ s 1 s 2 (s 1 - s 2 , 
C(s 2 ; s, sx)2s (s - s 2 ) 2 + si(s - 2s 1 + s 2 
+ C(si, s 2 , s)2s (s - si) 2 + s 2 (s + s x - 2s 2 

+ 5s 2 - 4s(si + s 2 ) - (si - s 2 ) 2 



, (D.2) 



The function £nc24 may be identified in the integrand on the right hand side of equation (2.11) 
in [31]. The quantity I4, which is defined there and contains higher order axial corrections to 



the Z width, is related to Q^c24 as follows: 

dsids 2 vA 



Sl s 2 3s ^(s- Sl ,s 2 ) = ---. 



(D.3) 



For numerical applications, the following limit is useful: 



lim^ NC24 (s; si, s 2 ) = 



9s 2 - A 2 



A 



A 2 3 (s 2 - A 



2^3 



(9s 4 + 6s 2 A 2 + A 4 ) + 0(A 2 )(I).4) 



with A = si — s 2 . 
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